X-Ray Resonant Scattering as a Direct Probe of Orbital Ordering in Transition-Metal 

Oxides 



On 
On 



O 

(N 
(N 



Michele Fabrizio (a) , Massimo Altarelli (6) , and Maurizio Benfatto (b ' c) 
(a) International School for Advanced Studies (SISSA) and INFM, Via Beirut 2-4, 34013 Trieste, Italy 
^ European Synchrotron Radiation Facility, Boite Postale 220, F-38043 Grenoble Cedex, France 
(c) INFN - Laboratori Nazionali di Frascati CP 13, 00044 Frascati, Italy 

(September 1997) 

X-ray resonant scattering at the K-edge of transition metal oxides is shown to measure the orbital 
order parameter, supposed to accompany magnetic ordering in some cases. Virtual transitions to the 
3d-orbitals are quadrupolar in general. In cases with no inversion symmetry, such as V2O3, treated 
in detail here, a dipole component enhances the resonance. Hence, we argue that the detailed 
structure of orbital order in V2O3 is experimentally accessible. 
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The orbital ordering in compounds containing ions 
with unfilled shells is a long standing problem |l| and 
has recently attracted new attention given the interest 
in the physics of transition metal oxides stimulated by 
cuprate superconductors and colossal magnetoresistance 
materials. In fact, while many experiments give direct 
information on magnetic structures, which are also con- 
sequence of partially empty shells, only indirect evidence 
points to the existence and structure of orbital order, un- 
less it is accompanied by a cooperative Jahn- Teller effect 
with a lattice distortion. An important and intriguing ex- 
ample is the prototype Mott-Hubbard insulator V2O3 [|| . 
At ambient pressure, this compound undergoes a metal- 
insulator transition at T c ~ 150K ||. The insulating 
phase below T c is antiferromagnetically ordered with a 
quite peculiar structure having both ferro- and antifcrro- 
magnetic bonds Q. Castellani, Natoli and Ranninger ||] 
argued that such a magnetic ordering implies an orbital 
ordering. Quite recently, Bao et al. Q have found by 
neutron diffraction that spin correlations in the metallic 
phase behave quite differently from the insulating one, 
giving further evidence of the importance of the orbital 
degrees of freedom. 

Actually, Castellani et al. || showed that (at least) 
two orbital structures, with different order parameters 
and wavevectors, are compatible with the observed mag- 
netic one. Later, based on the exchange constants fit- 
ted to neutron scattering data, one of the two orbital 
structures came to be favored j7],^],[| ■ However, there 
is still no direct observation of the orbital ordering, 
nor a firm experimental determination of its wavevec- 
tor, and of the temperature at which it appears, which 
may not coincide with the Neel temperature. In addi- 
tion to V2O3, there are many other examples where, in 
spite of the absence of direct evidences, the orbital de- 
grees of freedom are believed to play an important role, 
as the quite popular colossal magnetoresistance material 
Lai_ x (Ca,Sr) x Mn0 3 . 

In this Letter, we argue that it is indeed possible to 



have direct experimental access to the orbital ordering by 
X-ray resonant scattering, and that this process should 
be sufficiently intense to be readily observable with mod- 
ern synchrotron radiation sources. 

First of all we notice that, contrary to what we are 
going to argue for the resonant scattering, the non reso- 
nant X-ray diffraction has extremely little sensitivity to 
orbital order. We have estimated the cross section for 
orbital order scattering in this case by using hydrogenic 
3d-wavefunctions, and we find the corresponding scatter- 
ing power at the orbital order wavevector to be less then 
0.1 electrons per unit cell. 

Let us now show how an X-ray resonant elastic scatter- 
ing experiment can detect orbital ordering in transition 
metal compounds. The energy of the incoming beam is 
supposed to resonate with the K-edge of the transition 
element. Virtual transitions to the conduction band d- 
like components are quadrupolar, and the corresponding 
transition operator at a given atom in the lattice takes 
the form: 
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which also defines the dimensionless absorption opera- 
tor F = ^ a DaSa- In the above equation, d\ a is the 
operator which creates an electron of spin a in the i 
d-orbital, while s a creates a hole in the ls-shell, and 
C2 = [eui 2 T2^ds I c) \/ lnh I (15w), with u> the photon fre- 
quency, and n £i fc the density in the incoming beam of 
photons with polarization e and wavevector k. The radial 
matrix element is defined by r2.ds — J r2 drXd{ r ) r2 Xs{ r ), 
in terms of the radial wavefunctions Xs and Xd- 

If the atom is at position R, a phase factor e lk R weights 
the contribution of this single process. In obvious nota- 
tion, the absorption operator at this site is denoted by 
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Fr. The transition probability amplitude of an elastic 
scattering in which the above absorption process is fol- 
lowed by the emission of a photon with the same energy 
but polarization e' and wavevector k' is given by 
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where |<?) is the ground state, |n) an intermediate state 
and T/h the inverse lifetime of the core-hole due to other 
decay processes which might occur before the emission 
process takes place. 

At resonance becomes through ([!]) 
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In reality the resonant energy is not sharply defined, but 
the width of the d-like bands involved in the ordering is 
comparable to typical monochromator resolutions defin- 
ing the incoming energy and to the lifetime broadening 
r ~ leV. Therefore the resonant scattering via a K- 
edge quadrupole transition is a measure of the Fourier 
component at momentum Q — k — k' of an operator 



6 = J2D„(e',k')DUe,k), 
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which is a combination of products of d-orbitals annihi- 
lation and creation operators. The scattering geometry 
defines uniquely such a combination. The form of the 
scattering operator (Q) is similar to that of the orbital 
order parameter, generally written as A = J2i j a ij < ^i^ji 
where ay depend on the system under consideration (see 
below the specific case of V2O3). For a suitable choice 
of the scattering geometry, i.e. of wavevectors and po- 
larizations, the scattering operator (Q) includes a com- 
ponent equal to the orbital order parameter. This com- 
ponent will result in the only non vanishing amplitude 
for momentum transferred Q equal to the orbital order 
wavevector. The treatment sketched here can be applied 
without modification to cases like LaMnOs, in which the 
quadrupolar scattering channel is the only available one. 
In the case of V2O3, as we shall see, the absence of inver- 
sion symmetry at V-sites in the insulating phase allows 
dipole transitions as well. We will consider this more 
complicated but interesting case in detail. 

According to the general wisdom, the V ions in this 
compound have an oxidation state V 3+ , thus containing 
two d-electrons [Q. The five d-orbitals are split into a 
lower triplet of t2 9 orbitals and an higher doublet e g . It 
is common to use a reference frame in which the z-axis 
coincides with the c-axis of the non primitive hexago- 
nal cell, while the cc-axis is parallel to one of the bond 



connecting the V-atoms in the honeycomb lattice of the 
ab-plane. In this reference frame, the three t 2g orbitals 
become: d\ = d 3z 2_ r 2 ,and 
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The di orbital points towards the only nearest neighbor 
V along the c-axis, thus forming a strong covalent bond 
which is filled. Hence the d\ orbital, being inert, does not 
participate the orbital ordering, which only involves the 
d 2 and g? 3 orbitals. We notice that these orbitals, because 
of the crystal field in the corundum structure, acquire 
also a component of p-like symmetry, which we need to 
identify. A V-atom is surrounded by a distorted oxygen 
octahedron and by four vanadium second neighbors, pro- 
viding a non symmetric environment, as described in Ref. 
||. The crystal field potential on the central vanadium 
due to this environment may be parametrized by 



v cf = v 1 (r 3 ,_ 3 - 13,3) + 1^3,0. 
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As a result, the modified d-orbitals involved in the orbital 
ordering are, apart from the normalization, 

d 2 ^ d 2 + rj{up y + vp x ) ; d 3 -c d 3 + rj(up x ~ vp y ), (7) 

where r\u and r\v (u 2 + v 2 = 1) are proportional to Vi and 
V2, respectively, and inversely proportional to the energy 
difference between the vanadium 3d-orbitals and the p- 
orbitals involved. In order to obtain an estimate of this 
hybridization, an Xa calculation in the Z + 1 approxi- 
mation for a V06 cluster was performed, and provided 
77 ~ 0.2, and u ~ —0.37. A test of the adequacy of 
this approximation for a quantitative estimate is given 
by comparing the energy difference between the lowest 
absorption peak related to the d — p hybridized levels 
and the main absorption edge, which is 20.7eV in the 
cluster calculation, and 20.1eV experimentally ||. Also 
the relative intensities are in satisfactory agreement with 
experiment. 

According to Ref. two orbital orderings are com- 
patible with the observed magnetic structure. For the 
first one, the staggered orbital parameter is proportional 
to the operator 



case (A) : = d\d 3 + d\d 2 - 
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In the hexagonal reciprocal lattice, the Q-vector of this 
orbital ordering is predicted to be different from that 
of the magnetic ordering. In fact, while the latter has 
maximum diffraction at Q = (0.5,0.5,0) and zero at 
Q=(0.5,0.5,l), for the former the situation is reversed. 
This makes the search for the orbital ordering easier. In 
Cartesian coordinates, the wave vector of the orbital or- 
der (A) is 



M = 2tt| 0,-^=,- 



(9) 
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where a is the bond length in the honeycomb lattice of 
the basal plane, and c is the distance between two nearest 
neighbor V-atoms in the z-direction. The hexagonal cell 
that we use contains three layers, which is compatible 
with the notations used by Moon |jj and McWhan and 
Remeika ||. 

The other orbital order compatible with the magnetic 
structure has an order parameter proportional to 



case (B) : = d\d 2 — d\d 3 . 



(10) 



For this case the situation is even simpler, since the or- 
dering is no more staggered in the basal plane, while it 
remains staggered along the c-direction. However, since 
the two V-atoms on the basal plane belonging to the same 
hexagonal cell have opposite values of the order parame- 
ter, the smallest Q at which we expect a diffraction peak 
is 
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Let us now apply the general treatment described 
above to this specific example. However, since the con- 
duction band has both d and p-like character, we have 
to extend the general formalism to include, besides the 
quadrupole, also a dipole component in the transition 
amplitude. 

First of all, let us consider the quadrupole channel al- 
ready discussed. Since only the d 2 and d 3 combinations 
of Eq.(j^) are involved in the orbital ordering, it is more 
useful to pick up from (|l]) only the part involving these or- 
bitals. All the other allowed transitions are not expected 
to give a contribution to the elastic scattering at that 
particular Q. Hence, the relevant quadrupole absorption 
operator simplify to 



y=f ^ E [(- 4 a - (?. M 3 k) 
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where 1 / a/1 + rj 2 weighs the d-component in the conduc- 
tion band [see Eq. (|7|)1 , and the matrices M are 
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Next, let us study the dipole channel. Using the same 
notations as above, the dimensionless dipole absorption 
operator at resonance F\ = P<j s v i s defined through 
the dimensional one by 



j=x,y,z a 



(13) 



where Pj a creates a spin-cr electron in the j p- 
orbital, and C\ = eujri, ps y // 2Trh/ (3w), where t"i jPS = 
/ r 2 drXp{r)rXs(r), in analogy with the definition of r 2jC t s . 



The component of the dipole absorption operator Fi 
which is sensible to the orbital ordering can be easily 
found through Eq.(0), namely 
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where the factor in front of the sum is the weight of 
the p-component in the conduction band wave functions, 
and, for shortness, we have introduced the two vectors 
v 2 = (v, u, 0) and v 3 = (u, —v, 0). 

The whole transition operator is the sum of the 
quadrupole and dipole components. However Eqs.(^|), 
(0) and (Q) are still correct, provided we take 
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As a result, the scattering cross section will include, be- 
sides pure quadrupole and dipole contributions, an inter- 
ference term. To further proceed, we need to extract from 
the term proportional to the orbital order parameter, 
which can be in general written as (Tf\ + T 22 — iTfy) A a . 
Here a = A,B refers to the two possible order parame- 
ters Eq.® and Eq.@, and T£, T 2 Q 2 , T" 2 indicate the 
contributions of the pure dipole, the pure quadrupole and 
the interference channels, respectively. For these polar- 
ization and wavevector dependent quantities, we find: 



rpA 

± 22 



rpA 



rpA 
1 12 



rpi5 

J 22 — 



rpB 
1 ll 



rpB 

1 12 



w 2 (l + 


-if) 






c 2 




1 + V 2 


c 2 


[( 


V 


oCj 




1+T] 2 


2ujC 2 






( 


c 2 








-rj 2 ) 




v 2 


c\ 




1+T] 2 


CI 


[( 


V 


cCi 



M 2 /?) (?M 3 fc) + (2 <-> 3)] , (16) 



(17) 



(e*M 2 k>) (e • v 3 ) + {e*- v 2 ) (SM 3 k) 
3) 



1 + yy 2 2ujC 2 



e*M 2 k^j (eM 2 k) - (2 +■ 
'* ■ v 2 ) (e • v 2 ) - (2 «-> 3) 



3) 



(18) 
(19) 

(20) 



(2 



(e*M 2 k'^ (e-v 2 )+(e'' 
<-> 3) 



• v 2 



eM 2 k 



(21) 



From the cluster calculation, we obtain the radial matrix 
elements for dipole and quadrupole transitions, which al- 
low us to establish that the relative weight of the two 
terms in Eq.(|l5|) is given by r]Ci/C 2 ~ 7.4. Therefore 
we expect the pure dipole term alone to reproduce the 
total cross-section within ~ 13% accuracy. Hence, for 
the purpose of the present paper, it will be sufficient to 
calculate just this term. 

We are now in position to describe the behavior of 
the cross section for orbital order scattering as a func- 
tion of experimental configuration. Since the scattering 
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is elastic and the transferred momentum is fixed, also the 
scattering angle 27 is determined by sin7 = Qc/(2uj), 
being to the photon frequency, and Q the modulus of 
the wavevector for the orbital ordering, which will be ei- 
ther (||) or (11). Taking the K-edge of the vanadium 
fiLj = 5A65eV, and the lattice parameters in the insu- 
lating phase a — 2.88A and c = 2.70A, we find the two 
scattering angles — 15.49° and 73 = 17.34°. If we as- 
sume not to detect the polarization of the emitted light, 
only three degrees of freedom for the scattering geometry 
remain. The first is the angle <p of rotation of the scat- 
tering plane around the transferred momentum Q, with 
4> = corresponding to the component of the incoming 
wavevector k normal to Q lying in the afo-plane. The 
other two free scattering parameters are related to the 
polarization of the incoming beam, and reduce to one for 
a linearly polarized light. 

Let us now estimate the scattering cross section. The 
radial dipole matrix element as obtained from the clus- 
ter calculation is ri iPS = 780ro, where r$ is the classical 
electron radius. We also take T — 0.8eV 0. The re- 
sulting cross section is very large, of order 10 3 rQ, which 
is not surprising in view of our assumption of a single 
resonating level (similarly large cross sections for mag- 
netic X-ray scatterings were estimated by Hannon et al. 
JiH , for the 3d to 4/ resonance of a rare earth with one 
/-hole per atom). The dependence of the cross section 
on the angle </> is shown in Fig.l for the orbital order (A) 
of Eq.(|) and in Fig.2 for the case (B) of Eq.@. 




FIG. 1. Orbital order (A) scattering cross section 
in arbitrary units for the cr-polarization (solid line) and 
7r-polarization (dashed line). 

In conclusion it was shown that elastic resonant X-ray 
scattering in an appropriate experimental configuration 
provides a direct probe of the orbital order parameter. 
Numerical estimates show that, in general, the effect 
should be observable in transition metal oxides at the 
metal K-edge. This is even more so in the particularly 
interesting case of V2O3, where the absence of inversion 
symmetry in the insulating phase allows dipole transi- 
tions to relevant intermediate states. Interestingly, the 
same dipole matrix elements are expected to enhance the 
resonant magnetic scattering from the antiferromagnetic 



structure, allowing to monitor the interplay of orbital and 
magnetic order in the same experiment. 
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FIG. 2. Orbital order (B) scattering cross section 
in arbitrary units for the cr-polarization (solid line) and 
7r-polarization (dashed line). 

Our results imply that resonant X-ray experiments can 
determine the temperature dependence of the orbital or- 
der parameter, as well as its wavevector. These quantities 
have so far eluded any direct measurement. 

One of us (M.F.) thanks ESRF, where most of this 
work was performed, for its hospitality 
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